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It is shown that a multilayer comprised of alternating thin superconducting and insulating layers 
on a thick substrate can fully screen the applied magnetic field exceeding the superheating fields H s 
of both the superconducting layers and the substrate, the maximum Meissner field is achieved at an 
optimum multilayer thickness. For instance, a dirty layer of thickness ~ 0.1 |im at the Nb surface 
could increase H s ~ 240 mT of a clean Nb up to H s ~ 290 mT. Optimized multilayers of NbaSn, 
NbN, some of the iron pnictides, or alloyed Nb deposited onto the surface of the Nb resonator 
cavities could potentially double the rf breakdown field, pushing the peak accelerating electric fields 
above 100 MV/m while protecting the cavity from dendritic thermomagnetic avalanches caused by 
local penetration of vortices. 


I. INTRODUCTION 

The maximum magnetic held H which can be screened 
by a superconductor in the vortex-free Meissner state has 
attracted much attention, both from the fundamental 
and applied perspectives 1-16 . This problem is particu¬ 
larly important for the Nb resonator cavities 2 which have 
extremely high quality factors Q(2K) ~ 10 10 —10 11 up to 
the breakdown fields ~ 200 mT at which the screen¬ 
ing current density J approaches the depairing limit 1 , 
Jd — H c /Xq, where H c — 200 mT is the thermodynamic 
critical held of Nb and Ao — 40 nm is the London pen¬ 
etration depth. The lack of radiation losses and vortex 
dissipation in the Nb cavities (clean Nb has the highest 
lower critical held H c i ~ 180 mT among type-II super¬ 
conductors) enables one to probe the high-held nonlinear 
quasiparticle conductivity 17 which can be tuned by alloy¬ 
ing the surface with impurities 18-20 . 

The screening held limit of Nb could be exceeded by us¬ 
ing s-wave superconductors with higher H c and the crit¬ 
ical temperature T c , but such materials are prone to the 
dissipative penetration of vortices at H ~ H c \ < H^[ b . 
To address this problem it was proposed to coat the Nb 
cavities with multilayers of thin superconductors (S) with 
high H c > H^ b separated by dielectric (I) layers 11 (see 
Fig. 1). This approach is based on the lack of ther¬ 
modynamically stable parallel vortices in decoupled S 
screens of thickness d s < A, which also manifests itself 
in a strongly enhanced H c \ in thin films predicted by 
Abrikosov 5,6 and observed in Refs. 7 10. 

The maximum held H m screened by N superconduct¬ 
ing layers of thickness d = Nd s A is limited by 
the superheating held H s of S-layers 11 , for example, 
H s ~ 0.84ff c = 454 mT for Nb 3 Sn at T < T c . Here the 
Meissner screening currents at H = H s become unstable 
with respect to infinitesimal perturbations of electromag¬ 
netic held and the order parameter, while the magnetic 
barrier for penetration of vortices vanishes 12-14 . This pa¬ 
per addresses the limits to which the maximum screening 
held can be increased by S layers with given d s , A and 
H s deposited onto a thick Nb substrate with given Ao 
and H s q. It is shown that: 1. The maximum H m can 



FIG. 1: Distribution of the magnetic field h(x) in a multilayer 
of thickness d. The dashed lines show thin dielectric layers 
(■di d s ) separating the superconducting layers. 


be reached at an optimum multilayer thickness d m which 
depends on the materials parameters of S layers and the 
substrate; 2. The optimized S-I-S multilayer can screen 
the held exceeding both superheating helds, H s and H s0 ; 
3. S-I-S multilayers arrest thermomagnetic avalanches 
caused by local penetration of vortices at defects and do 
not let them develop into global hux jumps, which oth¬ 
erwise quench the cavity at helds well below H m . 


II. SUPERHEATING FIELD, DEPAIRING 
CURRENT DENSITY AND INSTABILITY OF 
MEISSNER STATE 

Distribution of a low-frequency (hco <C fceT c ) rf mag¬ 
netic held h(x) in a multilayer can be described by the 
London equations A 2 ft" = h\ at 0 < x < d and Aghf,' = (12 
at x > d, where the prime denotes differentiation with 
respect to x. Given that h{x) in a multilayer calculated 
numerically from the Eilenberger equations is close to the 
result of the London theory 13 , I hrst disregard the non¬ 
linear Meissner effect 21 which makes A dependent on J 
at J ~ Jd . Solutions of the London equations are (see 
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also Refs. 15,16): 


h\(x) = H[( 1 — c)e X ! x + ce x ^ A ], 0 < x < d, 
h 2 {x) = bHe {d - x)/Xo 1 x>d , 


^2 ’ 


b = 


(1 + k)g 


(1) 

( 2 ) 

( 3 ) 


Here h(x ) and the rf electric field E = —g,oOjX 2 h' are 
continuous at x = d, k = (A — Ao)/(A + Ao), and g = 
exp(d/A), both for a single S film and for a stack of S-I 
layers with di -C d s . 

For a S-I bilayer, the low-field surface resistance R s 
is determined by the total Joule rf power, H 2 R S / 2 = 
(<t/2) f 0 d E 2 ( x)dx + (ao/2) f^° E 2 (x)dx + qi , where cr and 
(To are the quasiparticle conductivities in S layers and 
the substrate, respectively, and qi accounts for dielectric 
losses. Using Eqs. (l)-(3), and R s = g,QU) 2 crX 3 /2 in the 
local dirty limit 3 yields 22 

it* = (k + S giy + - *dk g 2 /x\ 

R s og 2 { 1 + k) 2 4^d. i o; 3 e"A 2 Agff 2 

(k + g 2 ) 2 + (A + A 0 ) 2 (fc + 5 2 ) 2 ’ 


where i? s and i? s o include the residual resistances 3 , and 
e = e' — ie" is the complex dielectric permeability of a 
low-loss I layer with e" -C e'. 

For A > Ao, the constant c in Eqs. (l)-(3) is posi¬ 
tive so the current density J(0) = h[( 0) = (1 — 2c)H/\ 
at the surface is lowered by the counterflow induced 
by the substrate 16 . The Meissner state in S layers 
and the substrate is stable with respect to infinitesi¬ 
mal perturbations 12-14 if the current densities are smaller 
than the respective depairing limits, J(0) = h[( 0) < 
H s /X and J{d) = h' 2 {d) < H s0 /X 0 . These conditions 
define the field region of the vortex-free state: 


H .>*£ 


- kX 


Shown in Fig. 2 is the H — d diagram where the Meissner 
region is stable below both curves defined by Eqs. (5). If 
H s > H s qXq /A, thin S layers with d< C A are stable but do 
not screen the magnetic field which is thus limited by H s0 
of the substrate. As d increases, H{d) increases until d 
reaches the crossing point m at which the maximum field 
H m exceeds both H s and H s q . The latter results from the 
stabilizing effect of counterflow on J(0) = (1 — 2c)H/X, 
where c(d) given by Eq. (3) decreases as d increases. As a 
result, H(d) decreases with d at d > d m , approaching the 
superheating field of the S layer at d A. The optimal 
thickness d m is defined by the equalities in Eq. (5) which 
give a quadratic equation for g = exp(d m /A). Hence, 


d m = Aln(/x + vV 2 + k), (6) 

where g = H s A/(A + Xq)H s o. Substituting d m back to 
one of Eqs. (5) yields the maximum screening field 


H„ 





1/2 


( 7 ) 



FIG. 2: Magnetic phase diagram of a multilayer. Meissner 
state is stable below both solid lines calculated from Eqs. 
(5) for k = 1/2 and H a0 = 0.5 H s . Here H(d) > H s at 
d > d c = ln(p + y/g? — k), and the crossing point m defines 
the optimum thickness d m at which H(d) is maximum. 


At H = H m and d = d m , both current densities at x = 0 
and x = d reach the depairing limits for the respective 
materials, but H m in Eq. (7) exceeds both H s and H s0 
if A > Ao, because it is the current density J ~ Jd but 
not the magnetic held which makes the Meissner state 
unstable 1,12-14 . The Meissner region lying below both 
curves in Fig. 2 is similar to that was calculated nu¬ 
merically in the London model assuming that the energy 
barrier for the vortex parallel to the surface vanishes 15,16 . 
The London approach 5,6 used in Refs. 15,16 has the well- 
known limitations, particularly at the fields H « H s for 
which the Bean-Livingston barrier is only a few vortex 
core diameters away from the surface. As a result, the 
London model cannot be used for the calculations of the 
superheating held, so the calculations of Refs. 15,16 only 
give the right order of magnitude for H m but can hardly 
be used for quantitative evaluation of H m . 

Equations (5) obtained from the condition of instabil¬ 
ity of Meissner currents have no ambiguity resulting from 
the evaluation of the surface barrier in the London model 
with the rigid core cutoff 10,16 which fails at H « 1L S be¬ 
cause of strong deformation of the vortex core at the 
surface 23 . Calculations based on the Ginzburg-Landau 
(GL) theory have shown that the instability at H = H s 
is caused by lateral perturbations with the wavelength 
~ e 3/4 AV 4 decaying over the length ~ \fX^ perpendicu¬ 
lar to the surface, where £ is the coherence length 12 . At 
T ss T c , the theory gives H s = f(n)H c where f{n) de¬ 
creases as the GL parameter k = A/£ increases 12 , so that 
H s ft: 1.2 H c at k S3 1 and H s = (1 + 1/v / 2k)v / 5J? c /3 at 
k /S> 1. At T <C T c , only the limiting value /(«) —> 0.84 
for k /?> 1 has been calculated 13,14 . 

Taking A = 120 nm and H s = 0.84 H c ~ 454 mT in Eq. 
(7) for a dirty NbaSn (Ref. 26) on Nb with H s0 = 240 
mT and Ao = 40 nm, yields H m = 507 mT at d m = 
1.1A = 132 nm. Here the Meissner state persists up to 
the held ~ 12% higher than H s of NbsSn, consistent with 
the London model of Ref. 15. For NbN films 28 with £ = 9 
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nm and A = 180 nm, we obtain d m ~ 79 nm and H m — 
288 mT, while for pnictides 24 ’ 25 with A ~ 200 nm and 
H c ~ 1 T, Eqs. (6)-(7) give H s = 0MH s ~ 840 mT 
and H m ~ 872 mT at d m = 1.78A = 356 nm. The case 
of the two-band superconductor MgB 2 with H c ~ 230 
mT is more complicated as the rf dissipation can occur 
at a smaller field H d — at which the screening 

current decouples two weakly coupled a and n bands 2 7 . 
For the typical ratio of coherence lengths, /£„■ ~ 0.2 — 
0.3, the band decoupling field H d ~ 50 mT is consistent 
with the rf breakdown field of 42 mT at 4 K observed on 
MgB 2 /Al 2 0 3 bilayers on the Nb substrate 8 . 

The enhancement of H m does not necessarily require I 
layers because Eqs. (6)-(7) also describe a dirty layer at 
the surface where A is increased due to a shorter mean 
free path £. For instance, a dirty Nb layer with i ~ 2 
nm has A ~ Ao^o/^) 1 ^ 2 — 180 nm and £ = (^o) 1 ^ 2 — 9 
nm. Using H s w 0.84 H c for k = Xq/£ = 20, yields d m = 
0.44A = 79 nm and H m = 1AAH C = 288 mT, the same 
as for the above case of NbN and some 20% higher than 
H s o = 240 mT of pure Nb. In the limit of (Ao/A) 2 — > 0, 
Eq. (7) gives H m = y 7 H 2 + H 2 0 = 1.465ff c = 293 mT. 

To see how the nonlinear Meissner effect can change 
Eq. (7), H m is also calculated from the GL theory at 
k 1 for which the y component of the magnetic vector 
potential A(x) satisfies the equation 12-14 


A 2 A" - A + £ 2 A 3 = 0. (8) 

Consider a dirty layer with constant A and £ at 0 < x < d, 
a substrate with constant Ao and £ 0 , and the same T c 
in both regions. The first integrals of Eq. (8) are: 
A 2 A' 2 — A 2 + £ 2 A 4 /2 = Cat0<a;<c? and A \A! 2 — 
A 2 +£qA 4 /2 = 0 at x > d. Continuity of A(x) and A'(x) 
at x = d, the boundary condition A'(0) = —H m , and 
the GL pairbreaking conditions A(0) = </>o/ 27T\/3£ and 
A(d) = c/)o/2ttV3^o under which the Meissner superflow 
becomes unstable 12-14 yield (2Tr\H m /<f> 0 ) 2 —5/18£ 2 = C\ 
at x = 0, and [27r/i(d)A/</>o] 2 — 1/3 £q + £ 2 /18£q = C\ 
and [27r/i(d)Ao/^o] 2 = 5/18£q at x = d. Here Ci = 
(27 t/ <fio) 2 C , and <j >o is the magnetic flux quantum. Ex¬ 
cluding Ci and h(d) and solving for H m reproduces Eq. 
(7) in which H s —> H s , and 


Hi 


U-?- + £- 

V 5^ 0 2 5? 4 


Hi 


(9) 


Here H s = \/5H c /3, and both T c and H c = 

^o/2 3//2 7tA^ = ^o/2 3 ^ 2 7tAoCo are independent of the mean 
free path, according to the Anderson theorem 14 . Two 
last terms in the parenthesis of Eq. (9) resulting from 
the nonlinear Meissner effect give a negligible contribu¬ 
tion to H m for a dirty layer with £ < ^o- 


III. PENETRATION OF VORTICES AT 
DEFECTS AND THERMOMAGNETIC 
STABILITY OF MULTILAYERS. 


The maximum screening field H rn at which the Messier 
state becomes absolutely unstable with respect to in¬ 
finitesimal perturbations can hardly be reached under 
realistic operating conditions which require that the ac¬ 
celerating resonator cavities remain stable with respect to 
penetration of vortices, strong transient electromagnetic 
perturbations of charged beams, and local field enhance¬ 
ment at surface defects. In the multilayer approach 11 I 
layers are instrumental to assure the necessary stability 
margin with respect to local penetration of vortices at 
surface defects, which can otherwise trigger thermomag- 
netic flux jumps 29,30 particularly at low temperatures 
T<T C and extremely high screening currents J ~ J d at 
which the cavities operate. Misinterpretation of this issue 
has lead to speculations that neither I layers nor the en¬ 
hancement of B c i is necessary, so a few ym thick Nb 3 Sn 
film coating of the Nb cavities could just be protected 
against penetration of vortices by the Bean-Livingstron 
surface barrier 15 . This assumption contradicts a vast 
body of experimental data on magnetization of high-K 
type-II superconductors for which inevitable materials 
or topographic defects at the surface cause premature 
local penetration of vortices at H c i < H < H s , or even 
H < H c i due to grain boundary weak links 31 . 

The maximum screening field H m at which the Meiss¬ 
ner state is stable with respect to penetration of vortices 
can be evaluated from Eqs. (6)-(7) with effective H s and 
H s o depending on the operating conditions. For instance, 
in a multilayer with h(d) < H c io — 180 mT at H = H m , 
a vortex entering through a defect in the S layer cannot 
penetrate further into the bulk Nb. Let a surface defect 
cause local penetration of vortices as the current density, 
J(0) = H'( 0) = (3H S / A reaches a fraction < 1 of J d in 
S layer with d s < A and If J(0) is not too close to 

J dl the London theory shows 6,11 that the energy barrier 
U = eoln(l /f3) per unit length of a vortex in the S layer 
at J = fd J d coincides with the bulk Bean-Livingston sur¬ 
face barrier 5,6 at H = (3H s h » H c 3 and (£/d s ) 2 1, 

where eo = <^o/47rp 0 A 2 . The criterion J(0) < J d /2 as¬ 
sures a reasonable protection against penetration of vor¬ 
tices caused by low-angle grain boundaries in polycrys¬ 
talline Nb 3 Sn or pnictides 31 , or local field enhancement 
at typical topographic defects in the Nb cavities 2 . 

At H s —> H s /2 and H s0 — > 170 mT, Eqs. (6)-(7) give 
H m = 278 mT at d m = 0.8A = 96 nm for Nb 3 Sn. If Nb 
can withstand the field H s o —> 200 mT observed on the 
best cavities 2,3 , the maximum screening field could reach 
H m = 295 mT at d m = 0.67A = 80 nm. Increasing /? by 
the materials refinements of Nb 3 Sn could push the peak 
fields over 300 mT. Pnictides with H c ~ 0.9 T, such as 
Bao. 6 Ko. 4 Fe 2 As 2 (Ref. 24) could provide H m = 426 mT 
at d m = 1.21 A = 242 nm, /3 = 1/2, and H s o —> 200 mT. 

To inhibit dissipative penetration of vortices, S film 
can be subdivided by I layers into N thinner S layers 
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FIG. 3: Propagating vortex loop (a) turning into the vortex- 
antivortex pair (b) in S layer with a surface defect (black dot) 
which lowers the field threshold of vortex penetration. 

with d s = d m /N. At h{d) < H c 10 , even if a vortex 
penetrates at a defect in the first S layer, it could not 
propagate into the next S layer and further in the bulk 
Nb where it can cause a thermomagnetic avalanche 30 . As 
H(t ) reaches the critical value (3H S at a week spot, a vor¬ 
tex line cannot penetrate parallel to the surface but first 
originates at a defect as a small semi-loop as depicted in 
Fig. 3. The vortex semi-loop expands under the action of 
the perpendicular Meissner current until it hits the I layer 
where most part of the dissipative vortex core disappears 
in a loss-free flux channel connecting two short vortices of 
opposite polarity. Because of the magnetic flux compres¬ 
sion in the I layer and the substrate, the energy of a per¬ 
pendicular vortex ~ d s e o ln(L/£) diverges with the lateral 
film size L, while the energy of the vortex-antivortex pair 
~ d s io ln(zt/£) grows with the distance u, where eo ~ eo 
(Ref. 33). This vortex-antivortex pair expands during 
the positive rf cycle and contracts and annihilates as H (t) 
changes sign. The upper limit of the pair size u m can be 
estimated neglecting the long-range vortex-antivortex at¬ 
traction described by the last term in the dynamic equa¬ 
tion, rjii/2 = — eo /u, where = </g/27r£ 2 p n is the 

viscous drag coefficient, and p n is the normal state resis¬ 
tivity. Hence, u m ~ 2(f>oH/ Xpui ~ Pfpn/V2 -kkpqXu at 
H = /3H S and w = 2i tv, giving u m ~ 4 pm and the rf 
power 32 q ~ dJ 2 ^/^ ~ fi 2 f 2 d(j)oB c p n /np$X 2 ~ 2 pW 
for Nb 3 Sn at p n = 0.2 /Aim, X = 100 ran, /3 = 1/2, 
/ = 0.84, k = 20, v = 2 GHz and d s /X = 0.2. Taking 
into account attraction of antiparallel vortices in S layer 
reduces u m and q. Penetration of vortex semi-loops at 
a defect appears more realistic than the model of a long 
vortex parallel to the surface for which the ill-defined no¬ 
tion of H cl = 0 associated with the magnetic flux trapped 
in I layer 15 is irrelevant to the rf dissipation 29 . 

IV. DISCUSSION 

Localization of the rf power in a thin S layer inhibits 
expansion of multiple vortex loops in the bulk and stops 
dendritic thermomagnetic avalanches 30 that are partic¬ 
ularly pronounced at the extremely high screening cur¬ 


rent densities J ~ Jd and low temperatures T <C T c 
in the materials like Nb 3 Sn, NbN or pnictides with low 
p- 1 and thermal conductivity k. The multilayer thus 
significantly reduces vortex dissipation as compared to 
the bulk Nb 3 Sn, yet a thin Nb 3 Sn layer with d ~ 100 
ran may only slightly increase the thermal impedance 
of the cavity wall, G = a]} + d s /k s + di/ki + dub/k^b , 
where ax is the Kapitza thermal resistance. For dxb = 3 
mm, kxb — 20 W/mK, ax = 2 kW/m 2 K, the Nb 3 Sn 
multilayer with d s = 100 nm, k s ~ 10 -2 W/mK, and 
Al20 3 dielectric layers with di = 4 nm and ki = 0.3 
W/mK (Ref. 34) increases G by only ~ 5%. A thicker 
Nb 3 Sn him with d ~ 2 — 3pm doubles G and reduces 
the held of thermal quench 3 , in addition to the smaller 
H c i = (</>o/47rA 2 )(ln n + 0.5) < 130 rnT of bulk Nb 3 Sn 
with A > 65 nm (Ref. 26) as compared to H^[ b ~ 180 
mT. Theoretically, I layers provide the strongest pinning 
of propagating vortices and stop them more efficiently 
than holes in thin films which have been used to termi¬ 
nate dendritic hux avalanches 35 . 

The optimum number of S layers for particular mate¬ 
rials is determined by a balance between reduced vortex 
dissipation and suppression of superconductivity at the 
S-I interfaces. Here H s of ideal S layers with d s > (£A) 1//2 
remains the same as in the bulk 12 , contrary to the 
assertion 15 that H s is reduced at small d s . This claim 
was based on the artifacts of the London model discussed 
above and on taking into account only one right vortex 
image in Fig. 4 of Ref. 15 instead of summing up an 
inhnite chain of vortex-antivortex image dipoles which 
ensure that vortex currents do not cross the film sur¬ 
face. If this effect is properly taken into account 6,11 H s 
in a thin him (d < A) is the same as in a thick him 
(d > A). Moreover, the GL simulations of Ref. 15 for 
a Nb 3 Sn him on Nb show that H m (d) reaches the max¬ 
imum H m (d m ) ~ 1.08 H s at d m ~ 1.15A but remains 
larger than H s in the whole region A < d < 2X for which 
the numerical results were presented. 

High-held rf performance of the Nb cavities can be 
boosted by depositing not only materials with higher H s 
but also alloyed Nb-I-Nb multilayers which can increase 
H m and beneht from a signihcant raise of Q(H) with H 
in a wide held region 17-20 . A polycrystalline Nb mul¬ 
tilayer may be tuned by alloying and heat treatment to 
reduce the residual resistance 3,17,18 , and is also less prone 
to the current-blocking grain boundaries than the A-15 
or pnictide compounds 31 . Enhancement of the vortex 
penetration held by a dirty Nb/Al 2 0 3 bilayer deposited 
onto the Nb cavity was observed in Ref. 36. 

In conclusion, optimized multilayers can significantly 
increase the Meissner screening held while inhibiting dis¬ 
sipative penetration of vortices. Implementation of such 
multilayer coatings could potentially double the accel¬ 
erating held gradients of superconducting resonators as 
compared to the existent high-performance Nb cavities. 

This work was supported by DOE HEP under Grant 
No. DE-SC0010081. 
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